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We examine the properties of embedded clusters within 1 kiloparsec using new data from 
the Spitzer Space Telescope, as well as recent results from 2MASS and other ground-based 
near-infrared surveys. We use surveys of entire molecular clouds to understand the range and 
distribution of cluster membership, size and surface density. The Spitzer data demonstrate 
clearly that there is a continuum of star-forming environments, from relative isolation to dense 
clusters. The number of members of a cluster is correlated with the cluster radius, such that 
the average surface density of clusters having a few to a thousand members varies by a factor 
of only a few. The spatial distributions of 5p;/zer-identified young stellar objects frequently 
show elongation, low density halos, and sub-clustering. The spatial distributions of protostars 
resemble the distribution of dense molecular gas, suggesting that their morphologies result 
directly from the fragmentation of the natal gas. We also examine the effects of the cluster 
environments on star and planet formation. Although Far-UV and Extreme-UV radiation from 
massive stars can truncate disks in a few million years, fewer than half of the young stars in 
our sample (embedded clusters within 1 kpc) are found in regions of strong FUV and EUV 
fields. Typical volume densities and lifetimes of the observed clusters suggest that dynamical 
interactions are not an important mechanism for truncating disks on solar system size scales. 



1. INTRODUCTION 

Since PP IV, there have been significant advances in ob- 
servations of young stellar clusters from X-ray to millime- 
ter wavelengths. But while much of the recent work has 
concentrated on the stellar initial mass function (IMF) or 
protoplanetary disk evolution (e.g., Lada and Lada, 2003), 
less attention has been directed to discerning the structure 
of young embedded clusters, and the evolution of that struc- 
ture during the first few million years. Physical properties 
of young embedded clusters, such as their shapes, sizes, and 
densities, should inform theories of cluster formation. In 
this contribution, we describe recent results in which these 
properties are obtained for a representative sample of young 
(1-3 Myr), nearby (d < 1 kpc), embedded clusters. 

This contribution is motivated by three recent surveys 
made with the Spitzer Space Telescope: the Spitzer Young 
Stellar Cluster Survey - which includes Spitzer, near-IR, 
and millimeter-wave images of 30 clusters, the Spitzer 
Orion Molecular Cloud Survey - which covers 6.8 sq. 



degrees in Orion, and the Cores to Disks (c2d) Legacy 
program, which surveyed several nearby molecular clouds 
{Evans et ai, 2003). These surveys provide a comprehen- 
sive census of nearly all the known embedded clusters in 
the nearest kiloparsec, ranging from small groups of sev- 
eral stars to rich clusters with several hundred stars. A new 
archival survey from Chandra (ANCHORS) is providing X- 
ray data for many of the nearby clusters. Since PP IV, the 
Two Micron All Sky Survey (2MASS) has become widely 
used as an effective tool for mapping large regions of star 
formation, particularly in the nearby molecular clouds. This 
combination of X-ray, near-IR and mid-IR data is a power- 
ful means for studying embedded populations of pre-main 
sequence stars and protostars. 

Any study of embedded clusters requires some method 
of identifying cluster members, and we begin by briefly 
reviewing methods which have progressed rapidly since 
PP IV, including work from X-ray to submillimeter wave- 
lengths, but with an emphasis on the mid-infrared spectrum 
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covered by Spitzer. Beyond Section 2 we focus almost en- 
tirely on recent results from Spitzer, rather than a review 
of the literature. In Section 3, we discuss the cluster prop- 
erties derived from large-scale surveys of young embedded 
clusters in nearby molecular clouds, including their sizes, 
spatial distributions, surface densities, and morphologies. 
In Section 4 we consider the evolution of young embed- 
ded clusters as the surrounding molecular gas begins to dis- 
perse. In Section 5 we discuss theories of embedded cluster 
evolution, and in Section 6 consider the impact of the clus- 
ter environment on star and planet formation. Our conclu- 
sions are presented in Section 7. 

2. METHODS OF IDENTIFYING YOUNG STARS 
IN CLUSTERS 

2.1. Near- and Mid-infrared 

Young stellar objects (YSOs) can be identified and clas- 
sified on the basis of their mid-infrared properties {Adams 
et al, 1987; Wilking et al, 1989; Myers and Ladd, 1993). 
Here we review recent work on cluster identification and 
characterization based primarily on data from the Spitzer 
Space Telescope. 

Megeath et al. (2004) and Allen et al. (2004) devel- 
oped YSO classification schemes based on color-color di- 
agrams from observations taken with the Infrared Array 
Camera (IRAC) on Spitzer. Examining models of proto- 
stellar envelopes and circumstellar disks with a wide range 
of plausible parameters, they found that the two types of 
objects should occupy relatively distinct regions of the dia- 
gram. Almost all of the Class I (star+disk+envelope) mod- 
els exhibited the reddest colors, not surprisingly, with the 
envelope density and central source luminosity having the 
most significant effect on the range of colors. The Class II 
(star+disk) models included a treatment of the inner disk 
wall at the dust sublimation radius, which is a significant 
contributor to the flux in the IRAC bands. Models of the 
two classes generally occupy distinct regions in color space, 
indicating that they can be identified fairly accurately from 
IRAC data even in the absence of other information such as 
spectra. 

Comparison of these loci with YSOs of known types 
in the Taurus star forming region shows reasonably good 
agreement (Hartmann et al, 2005). Some degeneracy in the 
IRAC color space does exist; Class I sources with low enve- 
lope column densities, low mass infall rates or certain ori- 
entations may have the colors of Class II objects. The most 
significant source of degeneracy is from extreme reddening 
due to high extinction, which can cause Class II objects to 
appear as low-luminosity Class I objects when considering 
wavelengths A < 10 /j,m. 

The addition of data from the 24 /mi channel of the 
Multiband Imaging Photometer for Spitzer (MIPS) provides 
a longer wavelength baseline for classification, particularly 
useful for resolving reddening degeneracy between Class I 
and II. It is also crucial for robust identification of evolved 
disks, both "transition" and "debris", which lack excess 



emission at shorter wavelengths due to the absence of dust 
close to the star. Such 24 observations are limited, how- 
ever, by lower sensitivity and spatial resolution compared to 
IRAC, as well as the generally higher background emission 
seen in most embedded regions. Muzerolle et al. (2004) 
delineated Class I and II loci in an IRAC/MIPS color-color 
diagram of one young cluster based on the 3.6-24 spec- 
tral slope. 

The choice of classification method depends partly on 
the available data; not all sources are detected (or observed) 
in the 2MASS, IRAC, and MIPS bands. IRAC itself is sig- 
nificantly more sensitive at 3.5 and 4.5 /im than at 5.8 and 8 
/<m, so many sources may have IRAC detections in only the 
two shorter wavelengths, and require a detection in one or 
more near-IR bands to classify young stars (Gutermuth et 
al, 2004; Megeath et al, 2005; Allen et al, 2005). Guter- 
muth et al (2006) refined the IRAC+near-IR approach 
by correcting for the effects of extinction, estimated from 
the H — K color, and developed new classification criteria 
based on the extinction-corrected colors. 

It is useful to compare some of the different classifica- 
tion schemes. In Fig. 1 we plot first a comparison of the 
IRAC model colors from Allen et al (2004), Hartmann et 
al. (2005) and Whitney et al (2003). In general, the models 
predict a similar range of IRAC colors for both Class I and 
Class II sources. Also in Fig. 1 we plot the same sample 
of IRAC data (NGC 2068/71) from Muzerolle et al. (2006) 
in three color-color planes which correspond to the clas- 
sification methods discussed above. In all diagrams, only 
those sources with detections in the three 2MASS bands, 
the four IRAC bands, and the MIPS 24 fjm band were in- 
cluded. For the sake of comparison with pie-Spitzer work, 
the points are coded according to their K-24-fim SED slope. 
Prior to Spitzer, a commonly used 4-class system was de- 
termined by the 2-10 /im (or 2-20 /im) slope (a), in which 
a > 0.3 = Class I, -0.3 < a < 0.3 = "flat" spectrum, 
-1.6 < a < -0.3 = Class II, and a < -1.6 = Class III 
(photosphere) (Greene et al, 1994). A few of the sources 
in Fig. 1 have been observed spectroscopically and de- 
termined to be T-Tauri stars, background giants, or dwarfs 
unassociated with the cluster. These are indicated. The dia- 
grams also show the adopted regions of color space used to 
roughly distinguish between Class I and Class II objects. 

Classifications made with these methods are in general 
agreement with each other, though some differences are 
apparent. For example, roughly 30% of Class I objects 
identified with the Allen et al. method and detected at 24 
/zm appear as Class II objects in the IRAC/MIPS-24 color 
space, however many of these are borderline "flat spec- 
trum" sources where the separation between Class I and II is 
somewhat arbitrary and may not be physically meaningful. 

These classification methods implicitly assume that all 
objects that exhibit infrared excess are YSOs. However, 
there can be contamination from other sources, including 
evolved stars, AGN, quasars, and high-redshift dusty galax- 
ies. Since most of these unrelated objects are faint high- 
redshift AGN (Stern et al, 2005), we have found that a 
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Fig. 1. — Identifying and classifying young stars using near- and mid-infrared measurements. In the panel at top left, 
a comparison of predicted IRAC colors from Allen et al. (2004) (A04), Hartmann et al., (2005) (H05) and Whitney et 
al., (2003) (W03). Triangles represent Class II models with T c ff = 4000K and a range of accretion rates, grain size 
distributions, and inclinations. Squares and circles are Class 1/0 models for a range of envelope density, centrifugal radius, 
and central source luminosity. In the remaining panels, we plot the data for the embedded cluster NGC2068/71 (Muzerolle 
et al., 2006). Point types are coded according to the measured SED slope between 2 and 24 fim. Spectroscopically 
confirmed T-Tauri, giant, and dwarf stars are indicated. In the top right panel, the large rectangle marks the adopted domain 
of Class II sources; the Class I domain is above and to the right (adapted from Allen et ai, 2004). In the bottom right panel 
(Gutermuth et al, 2006), dereddened colors are separated into Class I and II domains by the dashed line. Diagonal lines 
outline the region where most of the classifiable sources are found. In the bottom left panel, the approximate domains of 
Class I and II sources are indicated by the solid lines. The dotted line represents the adopted threshold for excess emission 
at 3.6 and 4.5 /im; sources below this that exhibit large [8]-[24] excess are probably disks with large optically thin or 
evacuated holes (adapted from Muzerolle et ai, 2004). Arrows show extinction vectors for Ay = 30 (Flaherty et al., 
2006). These figures show that the various color planes considered here yield similar results when used to classify Spitzer 
sources. 
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magnitude cut of to 3 . 6 < 14 will remove all but approxi- 
mately 10 non-YSOs per square degree within each of the 
IRAC-only Class I and Class II loci, and all but a few 
non-YSOs per square degree from the IR AC/MIPS -24 loci, 
while retaining most if not all of the cluster population. 

2.2. Submillimeter and Millimeter 

The youngest sources in star forming regions are char- 
acterized by strong emission in the sub-millimeter and far- 
infrared, but ususally weak emission shortward of 24 /xm. 
These "Class 0" objects were first discovered in sub-mm 
surveys of molecular clouds (Andre et al, 1993). They are 
defined as protostars with half or more of their mass still 
in their envelopes, and emitting at least 0.5% of their lu- 
minosity at submillimeter wavelengths. Motivated in part 
by the discovery of Class objects, observers have im- 
aged many embedded clusters in their dust continuum emis- 
sion at millimeter and submillimeter wavelengths, reveal- 
ing complex filamentary structure and many previously un- 
known sources (e.g., Nutter et al, 2005; Sandell and Knee, 
2001; Motte et al, 2001, 1998). 

These submillimeter and millimeter wavelength images 
generally have tens to hundreds of local maxima, but only 
a small fraction of these are "protostars" having an inter- 
nal heating source; the rest are "starless cores" having a 
maximum of column density but no internal heating source. 
The standard way to determine whether a submm source is 
a protostar or a starless core is to search for coincidence 
with a infrared point source, such as a Spitzer source at 
24 or 70 /xm, or a radio continuum point source, such as 
a VLA source at 6 cm wavelength. For example the proto- 
stars NGC1333-IRAS 4A, 4B, and 4C in Fig. 2 are each 
detected at 850 /xm (Sandell and Knee, 2001), and each 
has a counterpart in VLA observations (Rodriguez et al., 
1999) and in 24 /im Spitzer observations, but not in the 
IRAC bands. In a few cases, Class protostars such as 
VLA 1623 have been identified from their submm emission 
and their radio continuum, but not from their mid-infrared 
emission, because their mid-infrared emission is too heavily 
extinguished (Andre et ai, 2000). 

2.3. X-ray 

Elevated X-ray emission is another signature of youth: 
young stellar objects have typical X-ray luminosity 1000 
times that of the Sun. The e-folding decay time for this X- 
ray luminosity is a few 100 million years (see e.g., Micela 
et ai, 1985; Walter and Barry, 1991; Dorren et ai, 1995; 
Feigelson and Montmerle, 1999). Although the X-ray data 
of young stellar clusters will be contaminated by AGN and 
other sources, this contamination can be reduced by iden- 
tifying optical/infrared counterparts to the X-ray sources. 
X-ray sources where the ratio of the X-ray luminosity to 
the bolometric luminosity (Lx/Lboi) ranges from 0.1% to 
0.01% are likely pre-main sequence stars. In contrast to the 
infrared techniques described in 2.1, which can only iden- 
tify Class 0/1 and II sources; X-ray observations can readily 
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Fig. 2.— IRAC 4.5 /xm and MIPS 24 /xm images of IRAS- 
4 in NGC 1333. MIPS detects each of the three VLA 
sources, while IRAC detects their outflows but not the driv- 
ing sources. 

detect class II and class III objects, with perhaps some bias 
toward class III objects (Flacomio et al, 2003). The main 
limitation of X-ray observations is the lack of sensitivity 
toward lower mass stars. A complete sample of stars re- 
quires a sensitivity toward souces with luminosities as low 
as 10 27 erg cm -2 s _1 (Feigelson et al, 2005), the sensi- 
tivity of most existing observations are an order of magni- 
tude higher. The observed X-ray luminosity is also affected 
by extinction. Depending on the energy of the source, the 
sensitivity can be reduced by a factor of ten for sources at 
A y <~10 (Wolket ai, 2006). 

2.4. Emission Lines and Variability 

Among other techniques for identifying young cluster 
members, spectroscopic surveys for emission lines and pho- 
tometric surveys for variability have been used success- 
fully at visible and near-IR wavelengths. The most com- 
mon means of identifying young stars spectroscopically is 
through detection of optical emission lines, in particular 
Ha at 6563 A (Herbig and Bell, 1988). Large-scale ob- 
jective prism (Wiramihardja et al, 1989; Wilking et al, 
1987) and later, wide-field multi-object spectroscopy (e.g., 
Hillenbrand et al, 1993) has been effective in identifying 
young stars in clusters and throughout molecular clouds, 
however they miss the deeply embedded members that are 
optically faint or invisible. This problem is partly alleviated 
by large-scale surveys for photometric variability in the op- 
tical and near-IR. Recent near-IR surveys by Kaas (1999) 
and Carpenter et al. (2001, 2002) have been successful at 
identifying young cluster members in Serpens, Orion and 
Chamaeleon, respectively. 

2.5. Star Counts 

Much of the work on the density, size, and structure 
of embedded clusters has relied on using star counts; in- 
deed, the distribution of 2.2 /xm sources were used to iden- 
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tify clusters in the Orion B cloud in the seminal work of 
Lada et al. (1991). Instead of identifying individual stars as 
members, methods based on star counts include all detected 
sources and employ a statistical approach toward member- 
ship, in which an average density of background stars is 
typically estimated and subtracted out. In this analysis, the 
star counts are typically smoothed to produce surface den- 
sity maps; a variety of smoothing algorithms are in the liter- 
ature (Gomez et al, 1993; Gladwin et al, 1999; Carpenter, 
2000; Gutermuth et al., 2005; Cambresy et al, 2006) 

The degree of contamination by foreground and back- 
ground stars is the most significant limitation for star count 
methods, and the efficacy of using star counts depends 
strongly on the surface density of contaminating stars. In 
many cases, the contamination can be minimized by setting 
a if -band brightness limit (Gutermuth, 2005; Lada et al, 
1991). To estimate the position dependent contamination 
by field stars, models or measurements of the field star den- 
sity can be combined with extinction maps of the molecular 
cloud (Carpenter, 2000; Gutermuth et al, 2005; Cambresy, 
2006). These maps are subtracted from the surface density 
of observed sources to produce maps of the distribution of 
embedded stars; however, these maps are still limited by the 
remaining Poisson noise from the subtracted stars. 

Star count methods have the advantage that they do not 
discriminate against sources without infrared excess, bright 
X-ray emission, variability, or some other indication of 
youth. On the other hand, they only work in regions where 
the surface density of member stars is higher than the sta- 
tistical noise from contaminating field stars. In Fig. 3 we 
show maps of the IRAS 20050 cluster derived from the 
if -band star counts and from the distribution of infrared- 
excess sources. In the case of IRAS 20050, we find that the 
star count method provides a better map of the densest re- 
gions (due in part to confusion with bright nebulosity and 
sources in the Spitzer data), while the lower density regions 
surrounding these peaks are seen only in the distribution of 
Spitzer identified infrared excess sources (due to the high 
density of background stars). 

3. THE STRUCTURE AND EVOLUTION OF CLUS- 
TERS: OBSERVATIONS 

3.1. Identifying Clusters in Large Scale Surveys of 
Molecular Clouds 

Unlike gravitationally bound open clusters or globular 
clusters, embedded clusters are not isolated objects. In most 
cases, molecular cloud complexes contain multiple embed- 
ded clusters as well as distributed populations of relatively 
isolated stars. Recent large scale surveys and all sky cat- 
alogs are now providing new opportunities to study the 
properties of embedded clusters through surveys of entire 
molecular clouds. The advantage of studying clusters by 
surveying entire molecular clouds is twofold. First, the sur- 
veys provide an unbiased sample of both the distributed and 
clustered populations within a molecular cloud. Second, the 
surveys result in an unbiased measurement of the distribu- 



tion of cluster properties within a single cloud or ensemble 
of clouds. For the remainder of this discussion, we will use 
the word "cluster" to denote embedded clusters of young 
stars. Most of these clusters will not form bound open clus- 
ters (Lada and Lada, 1995). 

We now concentrate on two recent surveys for young 
stars in relatively nearby (<1 kpc) molecular clouds. Car- 
penter (2000) used the 2MASS 2nd incremental point 
source catalog to study the distribution of young stars in 
the Orion A, Orion B, Perseus and Monoceros R2 clouds. 
Since the 2nd incremental release did not cover the entire 
sky, only parts of the Orion B and Perseus clouds were 
studied. More recently, Spitzer has surveyed a number of 
molecular clouds. We discuss here new results from the 
Spitzer Orion Molecular Cloud Survey (Megeath et al, 
2006) and the Cores to Disks (c2d) Legacy program sur- 
vey of the Ophiuchus Cloud (Allen et al, 2006). We use 
these data to study the distribution of the number of cluster 
members, the cluster radius, and the stellar density in this 
small sample of clouds. 

The advantage of using these two surveys is that they 
draw from different techniques to identify populations of 
young stellar objects. The analysis of the 2MASS data 
relies on star counting methods (Section 2.5), while the 
Spitzer analysis relies on identifying young stars with 
infrared-excesses from combined Spitzer and 2MASS pho- 
tometry (Section 2.1; Megeath et al. 2006). The 2M ASS 
analysis is limited by the systematic and random noise from 
the background star subtraction, making the identification 
of small groups and distributed stars subject to large uncer- 
tainties. The Spitzer analysis is limited to young stars with 
disks or envelopes. A significant number of young stars in 
embedded clusters do not show excesses; this fraction may 
range from 20% to as much as 50% for 1-3 Myr clusters 
(Haischet al, 2001). 

Carpenter (2000) identified stellar density peaks more 
than six times the RMS background noise, and defined a 
cluster as all stars in a closed 2a contour surrounding these 
peaks. Megeath et al. (2006) defined clusters as groups 
of 10 or more IR-excess sources in which each member is 
within a projected distance of 0.32 pc of another member 
(corresponding to a density of 10 stars pc -2 ). Only groups 
of ten or more neighbors are considered clusters. The clus- 
ters identified in the Spitzer survey are shown in Fig. 4. 

3.2. The Fraction of Stars in Large Clusters 

It is now generally accepted that most stars form in clus- 
ters (Lada and Lada, 1995), but quantitative estimates of 
the fraction of stars which form in large clusters, small clus- 
ters, groups and relative isolation are still uncertain. Porras 
et al. (2003) compiled a list of all known groups and clus- 
ters with more than 5 members within 1 kpc of the Sun, 
while Lada and Lada (2003) compiled the properties of a 
sample of 76 clusters with more than 36 members within 2 
kpc. Although these compilations are not complete, they 
probably give a representative sample of clusters in the 
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Fig. 3. — IRAS 20050 surface densities derived from the statistical technique applied to all stars (left), and from identifying 
the infrared-excess sources (right). In the left panel, all sources having K<16 are plotted as a function of their position. 
Contours show the surface density of K-band sources, starting at 1450 pc~ 2 (5a above median field star density) and 
increasing at intervals of 750 pc~ 2 . In the right panel, sources with infrared excess emission are plotted, and contours of 
their surface density are shown for 10, 60, 160, 360, 760, and 1560 pc~ 2 . The statistical technique (left) yields a higher 
peak surface density (^6000 pc~ 2 at the center) than the IR-excess technique (^3000 pc~ 2 ), but the latter is more sensitive 
to the spatially extended population of young stars. 



nearest 1-2 kpc. In the sample of Porras et al. (2003), 80% 
of the stars are found in clusters with N sta r > 100, and 
the more numerous groups and small clusters contain only 
a small fraction of the stars (also see Lada and Lada, 2003). 

In Fig. 5, we plot the fraction of members from the 
2MASS and Spitzer surveys as a function of the number of 
cluster size. Following the work of Porras et al. (2003), we 
divide the distribution into four sizes: N star > 100, 100 > 
N s tar > 30, 30 > N star > 10, and N star < 10. The 
main difference from the previous work is that we include 
a bin for N sta r < 10; these we refer to as the distributed 
population. All of the observed molecular clouds appear 
to contain a distributed population. Carpenter (2000) esti- 
mated that the fraction of stars in the distributed population 
were 0%, 20%, 27%, and 44% for the Orion B, Perseus, 
Orion A and Mon R2 cloud, respectively, although the es- 
timated fraction ranged from 0-66%, 13-41%, 0-61% and 
26-59%, depending on the assumptions made in the back- 
ground star subtraction. In the combined Spitzer survey 
sample, the fraction of distributed stars is 32-1 1%, 26-24%, 
and 25-21 % for the Ophiuchus, Orion A and Orion B clouds 
respectively. The uncertainty is due to contamination from 
AGN: the higher fraction assumes no contamination, the 
lower number assumes that the distributed population con- 
tains 10 AGN for every square degree of map size. The 



actual value will be in between those; extinction from the 
cloud will lower the density of AGN, and some of the con- 
taminating AGN will be found toward clusters. In total, 
these measurements suggest that typically 20-25% of the 
stars are in the distributed population. 

There are several caveats with this analysis. The first is 
the lack of completeness in the existing surveys. Carpen- 
ter (2000) considered the values of ./V star as lower limits 
due to incompleteness and due to the masking of parts of 
the clusters to avoid artifacts from bright sources. Com- 
pleteness is also an issue in the center of the Orion Nebula 
Cluster (ONC) for the Spitzer measurements. Also, we have 
not corrected the Spitzer data for the fraction of stars which 
do not show infrared excesses, the actual number of stars 
may be as much as a factor of two higher (Gutermuth et al., 
2004). 

Another uncertainty is in the definition of the clusters. 
The clusters identified by these two methods are not en- 
tirely consistent. For example, in Orion A there is an un- 
certainty in the boundaries of the ONC. There is a large 
halo of stars surrounding this cluster, and the fraction of 
young stars in large clusters is dependent on whether stars 
are grouped in the ONC, in nearby smaller groups, or the 
distributed population. Both the 2MASS and the Spitzer 
data lead to an expansive definition of this cluster, extending 
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Fig. 4. — The spatial distribution of all Spitzer identified infrared excess sources from the combined IRAC and 2MASS 
photometry of Orion A (left), Orion B (right) and Ophiuchus (bottom center). The contours outline the Bell Labs 13 CO 
maps for the Orion A and B clouds (Bally et ai, 1987; Miesch and Bally, 1994), and an Ay map of Ophiuchus (Huard, 
2006). The small grey dots show all the detections in the Spitzer 3.6 and 4.5 /im bands with magnitudes brighter than 15th 
and uncertainties less than 0.15. The large grey dots are the sources with infrared excesses. The black circles and triangles 
are sources found in clusters using the method described in Section 3.1; the two symbols are alternated so that neighboring 
clusters can be differentiated. Note that there are two clusters in the Orion A cloud which are below the lower boundary of 
the Bell Labs map. Each of the clouds has a significant distributed population of IR-excess sources. 
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Fig. 5. — The distribution of the fraction of stars in clus- 
ters taken from Carpenter (2000) (circles) and the Spitzer 
surveys of Orion and Ophiuchus (diamonds). The Spitzer 
surveys show a range, depending on whether corrections are 
made for background AGN. In both the 2MASS and Spitzer 
surveys, the distributed population (N star < 10) accounts 
for more than 20% of the total number of stars. 

beyond the Orion Nebula and incorporating the OMC2/3 
and NGC 1977 regions, as well the L1641 North group for 
the 2MASS analysis. The resulting cluster contains a signif- 
icant number of stars in a relatively low stellar density envi- 
ronment far from the O-stars exciting in the nebula, which 
differs significantly from the environment of the dense core 
of the cluster embedded in the Orion Nebula. The treatment 
of the ONC is critical to this analisys: 50% (for the 2MASS 
sample) to 76% (for the Spitzer sample) of the stars in large 
clusters (N star > 100) are found in the ONC. 

A final caveat is that these results apply to the current 
epoch of star formation in the nearest kiloparsec. While 
the largest cluster within 1 kpc is the ONC with 1000-2000 
members, a growing number of young super star clusters, 
which contain many thousands of stars, have been detected 
in our Galaxy. Super star clusters may bridge the gap be- 
tween embedded clusters in then nearest kiloparsec, and the 
progenitors of the globular clusters which formed earlier in 
our Galaxy's history. Thus, the distribution of cluster sizes 
we have derived may not be representative for other regions 
of the Galaxy, or early epochs in our Galaxy's evolution. 

3.3. The Surface Density of Stars in Embedded Clus- 
ters 

In a recent paper, Adams et al. (2006) found a corre- 
lation between the number of stars in a cluster and the ra- 
dius of the cluster, using the tabulated cluster properties in 
Lada and Lada (2003). They found that the correlation is 
even stronger if only the 2MASS identified clusters from 



10000 







1000 


A v = 1 / 


X / 


z" 

100 


/ X 

/ M 

/ °yJ 


/ A v =1 


10 







0.1 1.0 10.0 

Radius (pc) 



Fig. 6. — N s t a r vs. cluster radius for the 2MASS survey 
(crosses) of Carpenter (2000), the Spitzer Ophiuchus and 
Orion surveys (triangles) of Megeath et al. (2006) and Allen 
et al. (2006), and the Spitzer young stellar cluster survey 
(diamonds). Lines of constant column density are shown 
for a column density for Ay = 1 and Ay = 10. The aver- 
age surface density of cluster members varies by less than 
an order of magnitude. 

Carpenter (2000) were used, in which case the parameters 
were derived in a uniform manner. The same correlation 
is seen in a sample of clusters defined by Spitzer identi- 
fied IR-excess sources. This correlation is shown for the 
2MASS and Spitzer samples in Fig. 6. This relationship 
shows that while N sta r varies over 2 orders of magnitude 
and the cluster radius (R c i U ster) varies by almost 2 orders of 
magnitude, the average surface density of cluster members 
(A s tar I 't R 2 c iuster) varies by less than one order of magni- 
tude. The lower surface (A v = 1) envelope of this corre- 
lation may result in part from the methods used to identify 
clusters. In particular, for the many clusters surrounded by 
large, low surface density halos of stars, the measured ra- 
dius and density of these clusters depends on the threshold 
surface density or spatial separation used to distinguish the 
cluster stars from those in the halos. We can convert the sur- 
face densities of members into column densities of mass by 
assuming an average stellar mass of 0.5 M Q . Assuming a 
standard abundance of hydrogen, and the typical conversion 
from hydrogen column density to Ay, we plot lines of con- 
stant Ay in Fig. 6. In this figure the clusters are bracketed 
by lines equivalent to Ay ~ 1 and Ay ~ 10. Interestingly, 
this result is similar to one of Larson's laws for molecular 
clouds, that the average column density of gas in molecu- 
lar clouds is independent of cloud size and mass {Larson, 
1985; see also the chapter by Blitz et al). 
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3.4. The Spatial Structure of Embedded Clusters 

One of the major goals of the Spitzer young stellar clus- 
ter and Orion surveys is to systematically survey the range 
of cluster morphologies by identifying the young stellar ob- 
jects with disks and envelopes in these clusters. An initial 
result of this effort is displayed for ten clusters in Fig. 7, 
which shows the surface density of IR-excess sources. In 
this section, we give a brief overview of the common struc- 
tures found in embedded clusters, both in the literature and 
in the sample of clusters imaged with Spitzer. We also dis- 
cuss ISO and Spitzer observations of the youngest objects 
in these regions, the Class I and sources. 

Many of the clusters shown in Fig. 7 appear elongated; 
this had also been evident in some of the earlier studies 
of clusters (Carpenter et al, 1997; Hillenbrand and Hart- 
mann, 1998). To quantify this asymmetry, Gutermuth et al. 
(2005, 2006) compared the distribution of stars as a func- 
tion of position angle to Monte Carlo simulations of circu- 
larly symmetric clusters, and demonstrated that the elon- 
gation is statistically significant in three of the six clusters 
in their sample. The elongation appears to be a result of 
the primordial structure in the cloud; for the two elongated 
clusters which have 850 ^m dust continuum maps, the elon- 
gation of the cluster is aligned with filamentary structure 
seen in the parental molecular cloud. This suggests that 
the elongation results from the formation of the clusters in 
highly elongated, or filamentary, molecular clouds. 

Not all clusters are elongated. Gutermuth et al. (2005) 
found no significant elongation of the NGC 7129 cluster, 
a region which also showed a significantly lower mean and 
peak stellar surface density than the more elongated clusters 
in his sample. Since the cluster was also centered in a cavity 
in the molecular cloud (see Section 5); they proposed that 
the lack of elongation was due to the expansion of the clus- 
ter following the dissipation of the molecular gas. However 
not all circularly symmetric clusters are easily explained by 
expansion; Gutermuth et al. (2006) find two deeply embed- 
ded clusters with no significant elongation or clumps, but 
no sign of the gas dispersal evident in NGC 7129. These 
two clusters, Cepheus A and AFGL 490, show azimuthal 
symmetry, which may reflect the primordial structure of the 
cluster. 

Examination of Fig. 7 reveals another common structure: 
low density halos surrounding the dense centers, or cores, 
of the clusters. With the exception of AFGL 490 and per- 
haps Cepheus A, all of the clusters in Fig. 7 show cores and 
halos. The core-halo structure of clusters has been studied 
quantitatively through azimuthally smoothed radial density 
profiles (Muench et al. 2003). Although these density pro- 
files can be fit by power laws, King models, or exponen- 
tial functions (Hillenbrand and Hartmann, 1998; Lada and 
Lada, 1995; Horner et al, 1997; Gutermuth, 2005), the re- 
sulting fits and their physical implications can be mislead- 
ing. As pointed out by Hartmann (2004), azimuthally aver- 
aged density profiles can be significantly steepened by elon- 
gation (Hartmann (2004) argues this for molecular cores, 



but the same argument applies to clusters). A more sophis- 
ticated treatment is required to study the density profiles of 
elongated clusters. 

It has long been noted that young stellar clusters are 
sometimes composed of multiple sub-clusters (Lada et al, 
1996; Chen et al, 1997; Megeath et al., 1996; Allen et al., 
2002; Testi, 2002). Clusters with multiple density peaks or 
sub-clusters were classified as heirarchical clusters by Lada 
and Lada (2003). In some cases it is difficult to distinguish 
between two individual clusters and sub-clusters within a 
single cluster. An example are the NGC 2068 and NGC 
2071 clusters in the Orion B cloud (Fig. 4). These appear as 
two peaks in a more extended distribution of stars, although 
the cluster identification method described in Section 3.1 
separated the two peaks into two neighboring clusters. In 
the sample of Gutermuth et al. (2005, 2006), clumpy struc- 
ture was most apparent in the IRAS 20050 cluster (also see 
Chen et al, 1997). In this cluster, the sub-clusters are aso- 
ciated with distinct clumps in the 850 [im map of the asso- 
ciated molecular cloud. This suggests that like elongation, 
sub-clusters result from structures in the parental molecular 
cloud. 

3.5. The Distribution of Protostars 

If the observed morphologies of embedded clusters re- 
sult from the filamentary and clumpy nature of the parental 
molecular clouds, then the younger Class 0/1 objects, which 
have had the least time to move away from their star forma- 
tion sites, should show more pronounced structures than the 
older, pre-main sequence Class II and Class III stars. Lada 
et al. (2000) found a deeply embedded population of young 
stellar objects with large K — L colors toward the ONC; 
these protostar candidates showed a much more elongated 
and clumpy structure than the young pre-main sequence 
stars in the Orion Nebula. Using the methods described 
in Section 2.1, we have identified Class 0/1 and II objects 
in clusters using combined Spitzer and ground-based near- 
IR photometry. In Fig. 8, we plot the distribution of class 
0/1 and II sources for four clusters in our sample. In the 
LI 688 and IRAS 20050 clusters, the protostars fall prefer- 
entially in small sub-clusters, and are less widely distributed 
than the Class II objects. In the Serpens and GGD 12-15 
clusters, the protostars are organized into highly elongated 
distributions. An interesting example containing multiple 
elongated distributions of protostars is the "spokes" cluster 
of NGC 2264, which shows several linear chains of proto- 
stars extending from a bright infrared source (Teixeira et al., 
2006). These chains, which give the impression of spokes 
on a wheel, follow filamentary structures in the molecular 
cloud. These data support the view that the elongation and 
sub-clustering are indeed the result of the primordial distri- 
bution of the parental dense gas. It is less clear whether the 
observed halos result from dynamical evolution or originate 
in situ in less active regions of star formation surrounding 
the more active cluster cores. The current data suggest that 
the halos are at least in part primordial; class 0/1 objects are 
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Fig. 7. — The distribution of infrared excess sources in ten clusters surveyed with Spitzer. The contours are at 1, 10 and 
100 IR-excess sources pc~ 2 . These data clearly show that clusters are not circularly symmetric, but are often elongated. 
Some of the clusters, such as IRAS 20050, show distinct clumpy structure, although much of the small scale structure seen 
in the highest contours is due to statistical fluctuations in the smoothing scale. The three most circularly symmetric clusters 
are Cepehus A, AFGL 490 and NGC 7129; the irregular structure in these clusters is due in part to statistical fluctuations 
in regions of lower surface density. 



observed in the halos of many clusters (Gutermuth et al, 
2004; Megeath et al, 2004). 

The spacing of protostars is an important constraint on 
the physical mechanisms for fragmentation and possible 
subsequent interactions by protostars. Kaas et al. (2004) 
analyzed the spacing of Class I and II objects identified in 
ISO imaging of the Serpens cluster. They calculated the 
separations of pairs of Class I objects, and found that the 
distribution of these separations peaked at 0.12 pc. In com- 
parison, the distribution of separations for Class II objects 
show only a broad peak at 0.2 to 0.6 pc; this reflects the 
more spatially confined distributions of protostars discussed 
in the previous section. Teixeira et al. (2006) performed a 
similar analysis for the sample of protostars identified in 
the spokes cluster of NGC 2264. The distribution of near- 
est neighbor separations for this sample peaked at 0.085 pc; 
this spacing is similar to the Jeans length calculated from 
observations of the surrounding molecular gas. 

Although the observed typical spacing of protostars in 
Serpens and NGC 2264 apears to be ~ 0. 1 pc, as shown in 
Fig. 8, dense groups of protostars are observed in both these 



regions (and others) in which several Class 1/0 sources are 
found within a region 0. 1 pc in diameter. This is the dis- 
tance a protostar could move in 100,000 years (the nominal 
protostellar lifetime) at a velocity of 1 kms -1 . This sug- 
gests that if the velocity dispersion of protostars is compa- 
rable to the turbulent velocity dispersion observed in molec- 
ular clouds, interactions between protostars may occur, par- 
ticularly in dense groupings. On the other hand, obser- 
vations of some dense star-forming clumps show motions 
through their envelopes much less than 1 kms -1 (Walsh et 
al. 2004). The densest grouping of protostars so far iden- 
tified in the Spitzer survey is found in the spokes cluster. 
One of the protostars in the spokes has been resolved into 
a small system of 10 protostars by ground-based near-IR 
imaging and by Spitzer IRAC imaging. These protostars 
are found in a region 10,000 AU in diameter. It is not clear 
whether these objects are in a bound system, facilitating in- 
teractions as the sources orbit within the system, or whether 
the stars are drifting apart as the molecular gas binding the 
region is dispersed by the evident outflows (Young et al., 
2006). It should be noted that this group of 10 protostars 
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appears to be the only such system in the spokes cluster. 
Thus, although dense groups of protostars are present in star 
forming regions, they may not be common. 

4. GAS DISRUPTION AND THE LIFETIME OF 
EMBEDDED CLUSTERS 

In the current picture of cluster evolution, star forma- 
tion is terminated when the parental gas has dispersed. An 
understanding of the mechanisms and time scales for the 
disruption of the gas is necessary for understanding the du- 
ration of star formation in clusters, the lifetime and even- 
tual fate of the clusters, and the ultimate star formation ef- 
ficiency achieved in a molecular cloud. 

The most massive stars have a disproportionate effect 
on cluster evolution. Massive O stars can rapidly disrupt 
the parental molecular cloud through their ionizating radi- 
ation. The effect of the disruption is not immediate; once 
massive stars form in a molecular core, star formation may 
continue in the cluster while the massive star remains em- 
bedded in an ultracompact HII region. Examples of clusters 
in this state within 1 kpc of the Sun are the GGD 12-15 and 
Mon R2 clusters. The timescale for the disruption of the 
core is equivalent to the lifetime of the ultracompact HII 
region {Megeath et al, 2002); this lifetime is thought to be 
<~ 10, 000 years for the solar neighborhood (Casussus et al, 
2000; Comeron and Torra, 1996). 

In our sample of nearby embedded clusters, most sys- 
tems do not contain O-stars. However, a number of par- 
tially embedded clusters in the nearest 1 kpc show evidence 
for significant disruption by B type stars. Due to the par- 
tial disruption of the clouds, the clusters in these regions 
are found in cavities filled with emission from UV heated 
polycyclic aromatic hydrocarbons {Gutermuth et al, 2004). 
The time scale for the disruption by B stars can be esti- 
mated using measurements of the ages of the clusters. In 
our survey of nearby regions, we have three examples of re- 
gions with such cavities: NGC 7129 (earliest member B2), 
IC 348 (earliest member B5) and IC5146 (earliest member 
B0-1), with ages 2 Myr, 3 Myr and 1 Myr, respectively (Hil- 
lenbrand et al, 1992; Hillenbrand, 1995; Luhman et al, 
2003; Herbig andDahm, 2002). The presence of large, UV 
illuminated cavities in these regions suggest that the non- 
ionizing far-ultraviolet radiation (FUV) from B-stars may 
be effective at heating and evaporating molecular cloud sur- 
faces in cases where intense FUV radiation from O-stars is 
not present. For example, in the case of NGC 7129, Mor- 
ris et al. (2004) find that the temperature at the molecular 
cloud surface has been heated to 700 K by the FUV radia- 
tion. Future work is needed to determine if the high temper- 
atures created by the FUV radiation can lead to substantial 
evaporative flows. 

In regions without OB stars, however, some other mech- 
anism must operate. An example is IRAS 20050. Based 
on SCUBA maps, as well as the reddening of the members, 
Gutermuth et al. (2005) found that the cluster is partially 
offset from the associated molecular gas, suggesting that 



the gas had been partially dispersed by the young stars Al- 
though this region contains no OB stars, it displays multiple 
outflows (Chen et al. 1997). Another example may be the 
NGC 1333 cloud, where Quillen et al. (2005) found evi- 
dence of wind-blown cavities in the molecular gas. In these 
regions, outflows may be primarily responsible for dissi- 
pating the dense molecular gas (e.g., Matzner and McKee 
2000). 

It is important to note that star formation continues dur- 
ing the gas dissipation process. Even when the gas around 
the main cluster has been largely disrupted (such is the case 
in the ONC, IC 348 and NGC 7129), star formation contin- 
ues on the outskirts of the cluster in regions where the gas 
which has not been removed. Thus, the duration of star for- 
mation in these regions appears similar to the gas dispersal 
time of <~ 1 — 3 Myr. Older clusters have not been found 
partially embedded in their molecular gas (Leisawitz et al, 
1989). 

5. EARLY CLUSTER EVOLUTION 

Theories of cluster formation are reviewed elsewhere in 
these proceedings (see the chapters by Ballesteros-Paredes 
et al and Bonnell et al). Here we will discuss the dynam- 
ical evolution of young clusters during the first few million 
years. 

Although most stars seem to form within clusters of 
some type (see Section 4), only about ten percent of stars 
are born within star-forming units that are destined to be- 
come open clusters. As a result, for perhaps 90 percent of 
forming stars, the destruction of their birth aggregates is an 
important issue. Star formation in these systems is not 100 
percent efficient, so a great deal of cluster gas remains in 
the system. This gaseous component leaves the system in a 
relatively short time (a few Myr - see above) and its depar- 
ture acts to unbind the cluster. At the zeroth level of under- 
standing, if the star formation efficiency (SFE) is less than 
50% , then a substantial amount of unbinding occurs when 
gas is removed. However, this description is overly simple. 
The stars in the system will always have a distribution of 
velocities. When gas is removed, stars on the high velocity 
tail of the distribution will always leave the system (even 
for very high SFE) and those on the extreme low velocity 
tail will tend to stay (even for low SFE). The fraction of 
stars that remain bound after gas removal is thus a smooth 
function of star formation efficiency (several authors have 
tried to calculate the function: see Adams, 2000; Body and 
Kroupa, 2003a, 2003b; Lada et al, 1984). The exact form 
of the bound fraction, /&(e), which is a function of SFE, de- 
pends on many other cluster properties: gas removal rates, 
concentation of the cluster, total depth of the cluster poten- 
tial well, the distribution functions for the stellar velocities 
(radial vs isotropic), and the spatial profiles of the gaseous 
and stellar components (essentially, the SFE as a function 
of radial position). At the crudest level, the bound frac- 
tion function has the form /& m ^/e, but the aforementioned 
complications allow for a range of forms. 
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Fig. 8. — The spatial distributions of Spitzer identified class 1/0 (dark circles) and Class II (light circles) objects in four 
clusters: L1688 in Ophiuchus, Serpens, Mon R2 and IRAS 20050. The Class 1/0 sources are often distributed along 
filamentary structures, while the Class II sources are more widely distributed. Many small groups of protostars are dense 
enough that interactions bewteen individual objects may occur. 



The manner in which a cluster spreads out and dissolves 
after its gas is removed is another important problem. Af- 
ter gas removal, clusters are expected to retain some stars 
as described above, but such systems are relatively short- 
lived. For example, consider a cluster with N = 100 in its 
early embedded phase, before gas removal. After the gas 
leaves, typically one half to two thirds of the stars will be- 
come unbound along with the gas. The part of the cluster 
that remains bound will thus contain only N = 30 - 50 stars. 
Small groups with N < 36 have relaxation times that are 
shorter than their crossing times (Adams, 2000) and such 
small units will exhibit different dynamical behavior than 
their larger counterparts. In particular, such systems will re- 
lax quickly and will not remain visible as clusters for very 
long. 

As more data are taken, another mismatch between the- 
ory and observations seems to be emerging: The theoretical 
calculations described above start with an established clus- 
ter with a well-defined velocity distribution function, and 
then remove the gaseous component and follow the evo- 
lution. Given the constant column density relationship for 
clusters (section 3.4), that the velocity of the stars are viri- 



alized, and assuming that 30% of the cluster mass is in stars 
(see Lada and Lada, 2003), then the crossing time for the 
typical cluster in our sample is <~ 1 Myr (although it can 
be shorter in the dense centers of clusters). As a result, in 
rough terms, the gas removal time, the duration of star for- 
mation, and the crossing time are comparable. This implies 
that partially embedded clusters may not have enough time 
to form relaxed, virial clusters; this in turn may explain in 
part the range of morphologies discussed in Section 3. 

6. EFFECTS OF CLUSTERS ON STAR AND PLANET 
FORMATION 

The radiation fields produced by the cluster environment 
can have an important impact on stars and planets formed 
within. Both the extreme, ionizing UV (EUV) and the 
far-UV (FUV) radiation can drive disk evaporation (Shu 
et al., 1993; Johnstone et al, 1998; Storzer and Hollen- 
bach, 1999; Armitage, 2000). In the modest sized clus- 
ters of interest here (100-1000 stars), the mass loss driven 
by FUV radiation generally dominates (e.g., Adams et al, 
2004), although EUV photoevaporation can also be impor- 
tant (Armitage, 2000; Johnstone et al, 1998; Shu et al, 
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1993; Storzer and Hollenbach, 1999). For clusters with 
typical cluster membership e.g., with N sta r = 300 (Section 
3.1), the average solar system is exposed to a FUV flux of 
G w 1000 - 3000 (Adams et al, 2006), where G = 1 cor- 
responds to a flux of 1.6 xl0~ 3 erg cm~ 2 s _1 . FUV fluxes 
of this magnitude will evaporate a disk orbiting a solar type 
star down to a truncation radius of about 50 AU over a time 
scale of 4 Myr. As a result, planet forming disks are rela- 
tively immune in the regions thought to be relevant for mak- 
ing giant gaseous planets. Forming solar systems around 
smaller stars are more easily evaporated for two reasons. 
First, the central potential well is less deep, so the stellar 
gravity holds less tightly onto the disk gas, which is more 
easily evaporated. Second, we expect the disk mass to scale 
linearly with stellar mass so that disks around smaller stars 
have a smaller supply and can be evaporated more quickly. 
With these disadvantages, M stars with 0.25 M Q can be 
evaporated down to 10 AU in 4 Myr with an FUV radia- 
tion field of G = 3000. In larger clusters with more massive 
stars, Adams et al. (2004) find that regions with strong FUV 
and EUV can affect disks around solar mass stars on solar 
system size scales, truncating an initially 100 AU disk to a 
radius of 30 AU in 4 Myr. 

A full assesment of the importance of UV radiation on 
disks needs to be informed by the observed properties of 
clusters. What fraction of stars in the Spitzer and 2MASS 
samples are found in clusters with significant EUV radia- 
tion fields? We use the presence of an HII region as an in- 
dicator of a EUV field. In the Spitzer sample (the Orion A, 
Orion B and Ophiuchus clouds) the two clusters with HII re- 
gions contain 45% of the IR-excess sources. In the 2MASS 
sample (Orion A, Orion B, Perseus and Monocerous R2), 
55% of the young stars are found in the four clusters with 
HII regions. Thus, a significant fraction of stars is found 
in clusters with HII regions. However, in both the 2MASS 
and Spitzer samples most of the stars found in clusters with 
HII regions are found in the ONC. The ONC has a radius 
of 4 pc and many of the low mass stars in this cluster are 
more than a parsec away from the massive stars, which are 
concentrated in the center of the cluster. Thus, the fraction 
of stars exposed to a significant EUV field appears to be 
less than 50%. However, a more systematic determination 
of this fraction should be made as data become available. 

In addition to driving photoevaporation, EUV radiation 
(and X-rays) can help ionize the disk gas. This effect is po- 
tentially important. One of the most important mechanisms 
for producing disk viscosity is through magneto-rotational 
instability (MRI), and this instability depends on having 
a substantial ionization fraction in the disk. One problem 
with this idea is that the disk can become too cold and the 
ionization fraction can become too low to sustain the tur- 
bulence. If the background environment of the cluster pro- 
vides enough EUV radiation, then the cluster environment 
can be important for helping drive disk accretion. 

Clusters can also have an affect on the processes of star 
and planet formation through dynamical interactions. This 
raises a variant of the classic question of nature vs nurture: 



are the properties of the protostars and the emergent stars 
influenced by interactions or are they primarily the result 
of initial conditions in a relatively isolated collapse? The 
numerical simulations of cloud collapse and cluster forma- 
tion {Bate et al, 2003; Bonnell et al., 2003; Bonnell et al, 
2004) predict that interactions are important, with the indi- 
vidual protostars competively accreting gas from a common 
reservoir as they move through the cloud, and dynamical in- 
teractions between protostars resulting in ejections from the 
cloud. 

We assess the importance of interactions given our cur- 
rent understanding of cluster structure. The density of 
clusters, and of protostars in clusters, suggest that if stars 
move with velocities similar to the turbulent gas velocity 
(~ 1 kms -1 ), interactions can occur in the lifetime of a pro- 
tostar (100,000 yr). Gutermuth et al. (2005) found typical 
stellar densities of 10 4 stars pc~ 3 in the cores of two young 
clusters. If the velocity dispersion is 1 kms _1 , most pro- 
tostars will pass within 1000 AU - the size of a protostellar 
envelope - of another star or protostars within a protostellar 
lifetime. The observed spacing of Class 1/0 sources dis- 
cussed in Section 3.2 also suggests that interactions can oc- 
cur in some cases. At these distances protostars could com- 
pete for gas or interact through collisions of their envelopes. 
Interestingly, recent data suggest that, at least in some clus- 
ters, the observed pre-stellar clumps that make up the initial 
states for star formation are not moving dynamically, but 
rather have subvirial velocities (Walsh et al, 2004; Peretto 
et al, 2006). If these clusters are typical, then interactions 
between protostars in clusters would be minimal. 

Given the observed surface densities of clusters, is it pos- 
sible that a cluster could result from the collapse of individ- 
ual, non-interacting pre-stellar cores (i.e. nature over nur- 
ture)? If the starting density profile of an individual star for- 
mation event can be modeled as an isothermal sphere, then 
its radial size would be given by r = GM*/2a 2 s=a 0.03 
pc (where we use a typical stellar mass of M* = 0.5M Q 
and sound speed a = 0.2 kms ~ 1 ). A spherical volume of 
radius R = 1 pc can thus hold about 37,000 of these smaller 
spheres (in a close-packed configuration). Thus, we can 
conclude that there is no a priori geometrical requirement 
for the individual star forming units to interact. 

Once a star sheds or accretes its protostellar envelope, 
direct collisions are relatively rare because their cross sec- 
tions are small. Other interactions are much more likely to 
occur because they have larger cross sections. For example, 
the disks around newly formed stars can interact with each 
other or with passing binaries and be truncated (Kobayashi 
and Ida, 2001; Ostriker, 1994). In rough terms, these stud- 
ies indicate that a passing star can truncate a circumstellar 
disk down to a radius r<i that is one third of the impact pa- 
rameter. In addition, newly formed planetary systems can 
interact with each other, and with passing binary star sys- 
tems, and change the planetary orbits (Adams and Laugh- 
lin, 2001). In a similar vein, binaries and single stars can 
interact with each other, exchange partners, form new bi- 
naries, and/or ionize existing binaries (McMillan and Hut, 
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1996; Rasio etal, 1995). 

To affect a disk on a solar system (40 AU) scale requires 
a close approach at a distance of 100 AU or less. Gutermuth 
et al. (2005) estimated the rate of such approaches for the 
dense cores of clusters. They estimate that for the typical 
density of 10 4 stars per pc~ 3 , the interaction time is 10 7 
years, longer than the lifetime of the cluster. For N-body 
models of the modest sized clusters of interest here (100- 
1000 members), the typical star/disk system is expected to 
experience about one close encounter within 1000 AU over 
the next ~ 5 Myr while the cluster remains intact; close 
encounters within 100 AU are rare (e.g., Adams et al, 2006; 
Smith and Bonnell, 2001). Given that lifetime of the cluster 
is less than 5 Myr, these models again indicate a minimal 
effect on nascent solar systems. 

7. CONCLUSIONS 

1. The Distribution of Cluster Properties: Systematic sur- 
veys of giant molecular clouds from 2MASS and Spitzer, as 
well as targeted surveys of individual clusters, are providing 
the first measurements of the range and distribution of clus- 
ter properties in the nearest kiloparsec. Although most stars 
appear in groups or clusters, in many star-forming regions 
there is a significant distributed component. These results 
suggest that there is a continuum of star-forming environ- 
ments from relative isolation to dense clusters. Theories of 
star formation must take into account (and eventually ex- 
plain) this observed distribution. The 2MASS and Spitzer 
surveys also show a correlation between number of member 
stars and the radii of clusters, such that the average surface 
density of stars varies by a factor of only <~5. 

2. The Structure of Young Stellar Clusters: Common 
cluster morphologies include elongation, low density ha- 
los, and sub-clustering. The observed cluster and molec- 
ular gas morphologies are similar, especially when only the 
youngest Class 1/0 sources are considered. This similar- 
ity suggests that these morphologies (except possibly ha- 
los) result from the distribution of fragmentation sites in the 
parental cloud, and not the subsequent dynamical evolution 
of the cluster. Consequently, the surface densities and mor- 
phologies of clusters are important constraints on models of 
the birth of clusters. 

3. The Evolution of Clusters: The evolution of clusters 
is driven initially by the formation of stars, and then later 
on by the dissipation of gas. Gas dissipation appears to be 
driven by different processes in different regions, including 
photoevaporation by extreme-UV from O stars, photoevap- 
oration by far-UV radiation from B stars, and outflows from 
lower mass stars. Much of the gas appears to be dissipated 
in 3 Myr which is a few times the crossing time and the du- 
ration of star formation in these clusters. With these short 
timescales, clusters probably never reach dynamical equi- 
librium in the embedded phase. The survival of clusters as 
the gas is dispersed is primarily a function of the size of 



the cluster, the efficiency of star formation, and the rate at 
which the gas is dispersed. 

4. The Impact of Clustering on Star and Planet Formation: 
Far-UV and Extreme-UV radiation from massive stars can 
effectively truncate disks in a few million years. Extreme 
UV radiation is needed to affect disks around solar type 
stars on solar system scales (< 40 AU) in the lifetime of the 
cluster. Within our sample of molecular clouds, fewer than 
50% of the stars are found in regions with strong extreme 
UV-fields. The observed spacing of protostars suggest that 
dynamical interactions and competitive accretion may oc- 
cur in the denser regions of the observed clusters. How- 
ever, evidence of sub-virial velocities of pre-stellar conden- 
sations in at least one cluster hints that these interactions 
may not be important. Given the densities and lifetimes of 
the observed clusters, dynamical interactions do not appear 
to be an important mechanism for truncating disks on solar 
system size scales. 
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